home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2001 May / SGI IRIX Base Documentation 2001 May.iso / usr / share / catman / p_man / cat3 / complib / DSISL.z / DSISL
Encoding:
Text File  |  1998-10-30  |  1.5 KB  |  67 lines
  1.  
  2.  
  3.  
  4. DDDDSSSSIIIISSSSLLLL((((3333FFFF))))                                                            DDDDSSSSIIIISSSSLLLL((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DSISL   - DSISL solves the double precision symmetric system A * X = B
  10.      using the factors computed by DSIFA.
  11.  
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  14.       SUBROUTINE DSISL(A,LDA,N,KPVT,B)
  15.  
  16. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  17.      On Entry
  18.  
  19.      AAAA DOUBLE PRECISION(LDA,N)
  20.         the output from DSIFA.
  21.  
  22.      LLLLDDDDAAAA INTEGER
  23.         the leading dimension of the array  A .
  24.  
  25.      NNNN INTEGER
  26.         the order of the matrix  A .
  27.  
  28.      KKKKPPPPVVVVTTTT INTEGER(N)
  29.         the pivot vector from DSIFA.
  30.  
  31.      BBBB DOUBLE PRECISION(N)
  32.         the right hand side vector.  On Return
  33.  
  34.      BBBB the solution vector  X . Error Condition
  35.  
  36.      AAAA division by zero may occur if  DSICO  has set RCOND .EQ. 0.0 or  DSIFA
  37.      has set INFO .NE. 0  .  To compute  INVERSE(A) * C  where  C  is a matrix
  38.      with  P  columns
  39.         CALL DSIFA(A,LDA,N,KPVT,INFO)
  40.         IF (INFO .NE. 0) GO TO ...
  41.         DO 10 J = 1, P
  42.         CALL DSISL(A,LDA,N,KPVT,C(1,J)) 10 CONTINUE LINPACK.  This version
  43.      dated 08/14/78 .  James Bunch, Univ. Calif. San Diego, Argonne Nat. Lab.
  44.      Subroutines and Functions BLAS DAXPY,DDOT Fortran IABS
  45.  
  46.  
  47.  
  48.  
  49.  
  50.  
  51.  
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.